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<N ' 

As we search for an ever deeper understanding of the structure of hadronic matter one 
of the most fundamental questions is whether or not one can make a connection to the 
underlying theory of the strong interaction, QCD. We build on recent advances in the chiral 
extrapolation problem linking lattice QCD at relatively large "light quark" masses to the 
' physical world to estimate the scalar polarizability of the nucleon. The latter plays a key role 

in modern relativistic mean-field descriptions of nuclei and nuclear matter (such as QMC) 
and, in particular, leads to a very natural saturation mechanism. We demonstrate that 
the value of the scalar polarizability extracted from the lattice data is consistent with that 
■^j- | needed for a successful description of nuclei within the framework of QMC. In a very real 

sense this is the first hint of a direct connection between QCD and the properties of finite 
nuclei. 
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t^j- | §1. Introduction 

o _ 

Over the century since Rutherford made the first discovery concerning nuclear 
structure, understanding the mechanism for nuclear saturation has been a central 
concern. Within the traditional framework of many-body theory based on two-body 
potentials, strong short-distance repulsion was the key. The importance of relativity 
was made clear within quantum-hadrodynamics (QHD), with the scalar density of 
nucleons (which determines the attractive scalar, cr-meson, coupling) growing more 
^ . slowly than the density (which in turn determines the repulsive vector, w-meson, 

coupling). 

Since the discovery of QCD as the fundamental theory of the strong interac- 
tion, numerous attempts have been made to derive the NN force within quark mod- 
els. Fewer attempts have been made to understand nuclear structure at the quark 
level. 1 -'' 2 ' The quark-meson coupling (QMC) model 3 ) -6 ) stands between the tradi- 
tional meson-exchange picture and the hard core quark models. It is a mean field 
model in the sense of QHD but with the a and to mesons coupling to confined quarks, 
rather than to point-like nucleons. In the simple case of infinite nuclear matter, after 
a considerable amount of work, one finds that the only effect of the internal, quark 
structure of the nucleons is that the u-nucleon coupling becomes density dependent. 
That is, the nucleon effective mass takes the form: 

M* N = M N - g a (a)a (14) 

where g (7 (a) ~ g a + gaaO~, g a is the coupling at zero density and g aa is known as the 
scalar polarizability. 3 )' 6 )~ 8 ) 
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If the MIT bag model is used to describe the quark confinement (as in the 
original QMC model) the scalar polarizability is negative, that is it tends to oppose 
the applied scalar field. Indeed, one finds g aa = — 0.11g a R, with R the bag radius. A 
similar result is found in the NJL model, 7 ) once it is modified to simulate the effect 
of confinement 9 ) - a change which also resolves the long standing question of chiral 
collapse in the NJL model. 

The crucial point is now that a negative scalar polarizability (associated with 
a readjustment of the internal, quark structure of the nucleon in response to an 
applied scalar field) is all that is needed to lead to the saturation of the binding 
energy of nuclear matter as a function of density. In contrast with the Walecka 
model (QHD), 10 ) where the cr-nucleon coupling really is a constant, the mean scalar 
and vector fields at nuclear matter density are much smaller because one is not 
relying on a difference between ^n^n and ip^ N ipN arising from relativistic motion of 
the bound nucleons. 

While the situation just described is somewhat simplistic as a description of 
nuclear saturation, it is certainly sufficient. Moreover, it has recently been shown 
that a systematic expansion of the QMC model as a function of density allows a direct 
comparison with the widely used Skyrme effective forces. n ) In that framework one 
can also include anti-symmetrization at the nucleon level, leading to new values of 
the a and uj coupling constants. With these values the expansion of the effective 
model as an effective nuclear force yields agreement at the 10% level. This provides 
a remarkable indication of the role of confined quarks in even low energy nuclear 
physics. 

§2. The Status of Chiral Extrapolation for the Mass of the Nucleon 

The challenge of solving a strongly coupled field theory with spontaneous symme- 
try breaking is central to understanding the strong interaction within the Standard 
Model. It will be quite a few years before full QCD simulations can be performed at 
physical light quark masses - with computation time scaling like ~ l/m|j' 6 , 12 ) quark 
masses are typically restricted to above 50MeV for accurate calculations with im- 
proved quark and gluon actions. The exception to this is to use staggered fermions, 
which have shown some results at 25MeV, 13 ) but this approach has a number of 
technical problems, including multiple pions (extra tastes) - which complicate the 
chiral extrapolation problem which we explore next. 

Given that one cannot directly calculate hadron properties at realistic quark 
masses, if one wishes to make any connection with experimental data it is necessary 
to find an appropriate method to extrapolate the properties calculated at large light- 
quark mass to the physical value. This is the chiral extrapolation problem, which is 
complicated by spontaneous chiral symmetry breaking in QCD. As we shall explain, 
this problem has recently been solved by a careful reformulation of the effective field 
theory. 

The power of effective field theory for strongly interacting systems has been 
apparent for more than 50 years, beginning with Foldy's inclusion of the anomalous 
magnetic moment of the proton in a derivative expansion of the electromagnetic 
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coupling of photons to baryons. Of course, the fact that nucleon form factors can 
be approximated by a dipole with mass parameter M 2 = 0.71 GeV 2 means that the 
radius of convergence of the formal expansion is less than 0.84 GeV. In the last 
twenty years our understanding of the symmetries of QCD have led to a systematic 
expansion of hadron properties and scattering amplitudes about the limits m w = 
and q = 0, known as chiral perturbation theory. For the pseudoscalar mesons, which 
are rigorously massless in the limit where the u, d and s quark masses are zero (the 
chiral SU(3) limit), this approach has been especially successful. Although, even 
here, serious questions have been raised over whether the formal expansion is really 
convergent for masses as large as m s . 14 )~ 16 ) 

For our immediate interest in the scalar polar izability of the nucleon, we shall 
concentrate on understanding the dependence of Mn on light quark mass m q . Given 
the Gell-Mann-Oakes-Renner (GMOR) relation, m 2 oc m q , we prefer to use m 2 , 
rather than m q , as the measure of explicit chiral symmetry breaking in the expres- 
sions below. (We stress, however, that wherever the GMOR relation fails, e.g. for 
m 2 , > 0.8 GeV 2 or for — > in QQCD (where there are logarithmic chiral correc- 
tions to GMOR), one should revert to an expansion written explicitly in terms of 
m q .) The formal chiral expansion for Mjv, in terms of m w , about the SU(2) chiral 
limit in full QCD is: 

Mn = a + a 2 m 2 + a^m\ + a 6 m^. + . . . 

+a Nw (m n ) + <TAn(™>n) , (2-1) 

where ctb-k is the self-energy arising from a Bn loop (with B = N or A). These N 
and A loops generate the leading and next-to-leading non-analytic (LNA and NLNA) 
behaviour, respectively. Note that for the present we have omitted the tadpole term 
which also contributes at NLNA order. It adds complications which are unimportant 
for the present discussion - for details we refer to Sec. 2.3 and to Ref. 18 ) 

These loop integrals contain ultraviolet divergences which require some regu- 
larisation prescription. The standard approach is to use dimensional regularisation 
to evaluate the self-energy integrals. Under such a scheme the NNtt contribution 
simply becomes aN^m^) — > clna^-t an d the analytic terms, ao and c^w. 2 , undergo 
an infinite renormalisation. The A contribution produces a logarithm, so that the 
complete series expansion of the nucleon mass about m n = is: 

M N = c + c 2 m 2 + C4IJ1I + c L NA"i| + c NLNA ml log m n + . . . , (2-2) 

where the have been replaced by the renormalised (and finite) parameters q. 
The coefficients of the low-order, non-analytic contributions are known: 19 ) -21 ) 

CLNA = ~ 327/2^' CNLNA = 32T7 2 1^ ilz • (2 ' 3) 

Although strictly one should use values in the chiral limit, we take the experimental 
numbers with g A = 1.26, f w = 0.093 GeV, the N — A mass splitting, A = 0.292 GeV 
and the mass scale associated with the logarithms will be taken to be 1 GeV. 

We stress that whereas Eq. (2-2) was derived in the limit m n /A <C 1, at just 
twice the physical pion mass this ratio approaches unity. Mathematically the region 
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m n 2 (GeV 2 ) 

Fig. 1. Fits to lattice data with dipole regulator for the self-energy loop N — ► Ntt (neglecting the 
N — ► Ait term for simplicity). The dashed lines show the expansion up to successive orders in 
run, from m 2 up to m% - these alternate about the solid curve as we go through even and odd 
terms. 



f» A is dominated by a square root branch cut which starts at = A. Using 
dimensional regularisation this takes the form: 22 ) 



-3 32 2 1 
32tt/2 25 9a 2t7 



(24 3 - 3m* 4) log ^ \-2(A*-ml)i 



A + V^ 2 - rni 



(2-4) 



for m n < A, while for > A the second logarithm becomes an arctangent. 
Clearly, to access the higher quark masses in the chiral expansion, currently of most 
relevance to lattice QCD, one requires a more sophisticated expression than that 
given by Eq. (2-2). 

2.1. Issues of convergence 

Forgetting for a moment the issues surrounding the Air cut, we note that the 
formal expansion of the N — > iV7r — ► N self-energy integral, ctntt, has been shown 
to have poor convergence properties. Using a sharp, ultra-violet cut-off, Wright 
showed 23 -* that the series expansion, truncated at 0(m^), diverged for m w > 0.4 
GeV. This already indicated that the series expansion motivated by dimensional 
regularisation would have a slow rate of convergence. A similar conclusion had been 
reached somewhat earlier by Stuckey and Birse. 17 ) 

In considering the convergence of the truncated series, Eq. (2-2), it is helpful 
to return to the general form from which it was derived, namely Eq. (2-1). The 
dimensionally regulated approach requires that the pion mass should remain much 
lighter than every other mass scale involved in the problem. This requires that 
"t-tt/^xSB *C 1 (and m^/A <C 1 if we use the simple logarithm in Eq.(2-2) rather 
than the full cut structure of Eq.(2-4) ). 

An additional mass scale, which we address in detail below, is set by the physical 
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extent of the source of the pion field. 24 ) This scale, which is of order Rsource^ 
corresponds to the transition between the rapid, non-linear variation required by 
chiral symmetry and the smooth, constituent- quark like mass behaviour observed in 
lattice simulations at larger quark mass. An alternative to dimensional regularisation 
is to regulate Eq. (2-1) with a finite ultra-violet cut-off in momentum space. 

General considerations in effective field theory, as discussed for example by LeP- 
age, 25 ) suggest that one should not take the cut-off to oo. That is, one should not 
use dimensional regularization. Nor should one attempt to determine the "true" 
cut-off of the theory. Rather, one should choose a cut-off scale somewhat below the 
place where the effective theory omits essential physics and use data to constrain 
the renormalization constants of the theory and hence eliminate the dependence on 
that cut-off as far as possible. 

The issue is then what mass scale determines the upper limit beyond which the 
effective field theory is applicable. This is commonly taken to be A x pt ~ 47r/ 7r ~ 1 
GeV. Unfortunately this is incorrect for baryons. In the context of nuclear physics 
it has long been appreciated that nuclear sizes could never be derived from naive 
dispersion relation considerations of nearest t-channel poles - anomalous thresholds 
associated with the internal structure of nuclei dominate. Similarly, the size of a 
baryon is determined by nonperturbative QCD beyond the scope of chiral pertur- 
bation theory. The natural scale associated with the size of a nucleon is the inverse 
of its radius or a mass scale Rsource ~ 0.2 to 0.5 GeV - far below A x pt- In the 
context of effective field theory it is inconsistent to keep loop contributions from 
momenta above this scale! 

As a result of these considerations, we are led to regulate the radiative corrections 
which give rise to the leading and next-to-leading non-analytic contributions to the 
mass of the nucleon using a finite range regulator with a mass in the range Rsource- 
As far as possible residual dependence on the specific choice of mass scale (and form 
of regulator function) will be eliminated by fitting the renormalization constants to 
nonperturbative QCD - in this case data obtained from lattice QCD simulations. 
The quantitative success of applying the method is to be judged by the extent to 
which physical results extracted from the analysis are indeed independent of the 
regulator (over a physically sensible range) . We refer to this approach as finite range 
regularization, FRR. 

2.2. Analysis of lattice QCD data 

In our analysis we have taken as input both the physical nucleon mass and recent 
lattice QCD results of the CP-PACS Collaboration 26 ) and the JLQCD Collabora- 
tion. 27 ) This enables us to constrain an expression for Mjy as a function of the quark 
mass. The lattice results of Ref. 26 ) have been obtained using improved gluon and 
quark actions on fine, large volume lattices with high statistics. These simulations 
were performed using an Iwasaki gluon action 28 ) and the mean-field improved clover 
fermion action. In this work we concentrate on only those results with m sca = m va \ 
and the two largest values of (i.e., the finest lattice spacings a ~ 0.09 - 0.13 fm). 
We use just the largest volume results of Ref., 27 ) where simulations were performed 
with non-perturbatively improved Wilson quark action and plaquette gauge action. 
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The lattice volumes and spacings are similar for the two data sets. 

We chose to set the physical scale at each quark mass using the UKQCD 
method. 29 ^ That is, we use the Sommer scale r$ = 0.5 fm. 30 )' 31 ) This choice is 
ideal in the present context because the static quark potential is insensitive to chiral 
physics. This ensures that the results obtained represent accurate estimates of the 
continuum, infinite-volume theory at the simulated quark masses. The lattice data 
lies in the intermediate mass region, with to 2 between 0.3 and 1.0 GeV 2 . 

Based on Eq. (2-1) we evaluate the loop integrals a^ n = CLNAlwimn, A) and 
o-Aw = c N LNAlnA{m n ,A), where: 



-I 

7T Jo 



and 



= r dk k ' u2{k \ (2. 6 ) 

3 Jo Vk 2 + m 2 (A + Vfc 2 + m 2 ) ' V ' 

We use four different functional forms for the finite-ranged, ultra-violet vertex regula- 
tor, u(k) — namely the sharp-cutoff (SC), O(A-k); monopole (MON), A 2 /(A 2 + k 2 ); 
dipole (DIP), A 4 /(A 2 + k 2 ) 2 ; and Gaussian (GAU), ex.p(-k 2 / A 2 ). Closed expres- 
sions for these integrals in the first three cases are given in the Appendix of Ref. 32 ) 
Provided one regulates the effective field theory below the point where new short 
distance physics becomes important, the essential results will not depend on the 
cut-off scale. 25 ** We use knowledge learned in Ref. 32 ) as a guide for the appropriate 
scales for each of the regulator forms. 

We have fit the set of lattice data discussed above with the chiral expansions 
based on six different regularisation prescriptions. In all schemes we allowed the 
coefficients of the analytic terms up to ra\ to be determined by the lattice data. 
The regulator masses are once again fixed to their preferred values. The best fits, 
shown in Fig. 2, serve to highlight the remarkable agreement between all finite- 
range parameterisations. This demonstrates that the extrapolation can be reliably 
performed using FRR with negligible dependence on the choice of regulator. 

In view of the general discussion presented earlier, it is important to check the 
residual dependence on the choice of cutoff scale A. The resulting variation of the 
extrapolated nucleon mass, at the physical pion mass, for dipole masses ranging from 
0.6 to 1.0 GeV is shown in Fig. 3. We see that the residual uncertainty introduced 
by the cutoff scale is less than 2%, which is insignificant compared to the statistical 
error in extrapolating such a large distance. With the present data this statistical 
error is found to be around 13%. 

2.3. Complete analysis of the nucleon mass 

Here we summarise the results of the most recent and complete analysis of Ref. 18 ) 
Including the tadpole term, the nucleon mass has the formal expansion: 

M N = a + a 2 m 2 + a 4 m 4 + a§m% + a% N + a^ NA + of ad . (2-7) 
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m n 2 (GeV 2 ) 



Fig. 2. Extrapolation of lattice data for various regularisations prescriptions, without the tadpole 
term. The four indistinguishable, solid curves show the fits based on four different FRRs. 
The short dash-dot curve corresponds to the dimensionally regulated case and the long dash-dot 
curve to that where the branch point in Eq. (2-4) is included in dimensional regularisation. 

1.1 , , 
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A (GeV) 

Fig. 3. Illustration of the relatively weak dependence of the extrapolated nucleon mass on the choice 
of dipole regulator mass, A. 



As explained earlier, within dimensional regularisation this leads to a truncated 
power series for the chiral expansion, 

M N = c + c 2 ml + clna"4 + c±m\ + c N lna"4 In — + c^m% + ... , (2-8) 

A* 

where the bare parameters, aj, have been replaced by the finite, renormalised coef- 
ficients, q. Through the chiral logarithm one has an additional mass scale, fi, but 
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the dependence on this is eliminated by matching C4 to "data" (in this case lattice 
QCD). We work to fourth order in the chiral expansion and include the next analytic 
term to compensate short distance physics contained in the NLNA loop integrals as 
suggested in Ref. 33 ) 

On the other hand, the systematic FRR expansion of the nucleon mass is: 

M N = ciq + a^ml + aiml + a§m% + a^ N (m w , A) 

+a n NA (m n ,A)+a^ d (m n ,A), (2-9) 

where the dependence on the shape of the regulator is implicit. The dependence 
on the value of A and the choice of regulator is eliminated, to the order of the 
series expansion, by fitting the coefficients, a^, to lattice QCD data. The clear 
indication of success in eliminating model dependence and hence having found a 
suitable regularisation method, is that the higher order coefficients (af , i > 4) should 
be small and that the renormalised coefficients, q, and the result of the extrapolation 
should be insensitive to the choice of ultraviolet regulator. The forms of the 7ViV7r 
and NAir contributions were given earlier, while the tadpole term has the form: 




(2-10) 



where in Eq. (2T0) to, defined such that the term in braces vanishes at m w = 0, is a 
local counter term introduced to ensure a linear relation for the renormalisation of 

The result of the fit to the latest two-flavor lattice QCD data for the mass of the 
nucleon is shown in Fig. 4, with the corresponding parameters given in Table I. It 
is remarkable that all four curves based on FRR are indistinguishable on this plot. 
Furthermore, we see from Table I that the coefficient of m\ in all of those cases 
is quite small - an order of magnitude smaller than the dimensionally regularised 
forms. Similarly, the FRR coefficients of m% are again much smaller than their DR 
counterparts. This indicates that the residual series, involving a^, is converging when 
the chiral loops are evaluated with a FRR. 



Regulator 


ao 


CS2 


04 


0,6 


A 


X 2 /dof 


Dim. Reg. 


0.827 


3.58 


3.63 


-0.711 




0.43 


Dim. Reg. (BP) 


0.792 


4.15 


8.92 


0.384 




0.41 


Sharp Cutoff 


1.06 


1.47 


-0.554 


0.116 


0.4 


0.40 


Monopole 


1.74 


1.64 


-0.485 


0.085 


0.5 


0.40 


Dipole 


1.30 


1.54 


-0.492 


0.089 


0.8 


0.40 


Gaussian 


1.17 


1.48 


-0.504 


0.095 


0.6 


0.40 



Table I. Bare, unrenormalised, parameters extracted from the fits to lattice data displayed in Fig. 4 
- from Ref. 18 ' All quantities are in units of appropriate powers of GeV and fi = 1 GeV in 
Eq. (2-8). 

The corresponding, physically meaningful, renormalized coefficients are shown in 
Table II, in comparison with the corresponding DR coefficients found using Eq. (2-8). 
Details of this renormalisation procedure are given in Ref. 32 ) Thede gree of consis- 
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0.0 0.2 0.4 0.6 0.8 1.0 

m n 2 (GeV 2 ) 

Fig. 4. Fits to lattice data for various ultra-violet regulators including the tadpole term - from 
Ref. 18 ' The sharp cut-off, monopole, dipole and Gaussian cases are depicted by solid lines, 
indistinguishable on this plot. The dimensional regularised forms are illustrated by the dash- 
dot curves, with the correct branch point corresponding to the higher curve. Lattice data is 
from Ref. 26) 



Regulator 


CO 


C2 


C4 


Dim. Reg. 


0.827(120) 


3.58(50) 


3.6(15) 


Dim. Reg. (BP) 


0.875(120) 


3.14(50) 


7.2(15) 


Sharp cutoff 


0.923(130) 


2.61(66) 


15.3(16) 


Monopole 


0.923(130) 


2.45(67) 


20.5(30) 


Dipole 


0.922(130) 


2.49(67) 


18.9(29) 


Gaussian 


0.923(130) 


2.48(67) 


18.3(29) 



Table II. Renormalised expansion coefficients in the chiral limit obtained from various regulator 
fits to lattice data - from Ref. 18 ^ (All quantities are in units of appropriate powers of GeV.) 
Errors are statistical in origin arising from lattice data. Deviations in the central values indicate 
systematic errors associated with the chiral extrapolation. 



tency between the best-fit values found using all choices of FRR is remarkably good. 
On the other hand, DR significantly underestimates C4. We can understand the 
problem very simply; it is not possible to accurately reproduce the necessary 
behaviour of the chiral loops (for > A) with a 3rd order polynomial in m\. 

It is clear that the use of an EFT with a FRR enables one to make an accurate 
extrapolation of the nucleon mass as a function of the quark mass. Indeed, all of the 
FRR used yield physical nucleon masses that lie within a range of 1%. It is especially 
interesting to observe the very small difference between the physical nucleon masses 
obtained with each FRR when we go from LNA to NLNA - i.e., when the effect of 
the A is included. (The change is typically a few MeV for a FRR but more than 
100 MeV for DR.) Once again the convergence properties of the FRR expansion are 
remarkable. 
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LNA 




NLNA 






Regulator 


m N 


mjv 


(J AT 


Dim. Reg. 


0.784 


0.884 


±0.103 


50.3 


± 


10.0 


Dim. Reg. (BP) 


0.784 


0.923 


±0.103 


42.7 


± 


10.0 


Sharp cutoff 


0.968 


0.961 


±0.116 


34.0 


± 


13.0 


Monopole 


0.964 


0.960 


±0.116 


33.0 


± 


13.0 


Dipole 


0.963 


0.959 


±0.116 


33.3 


± 


13.0 


Gaussian 


0.966 


0.960 


±0.116 


33.2 


± 


13.0 



Table III. The nucleon mass, tum (GeV), and the sigma commutator, on (MeV), extrapolated 
to the physical pion mass obtained in a NLNA (fourth order) chiral expansion - from Ref. 18 ' 
Convergence of the expansion is indicated by the nucleon mass obtained in an analysis where 
we retain only the LNA (third order) behaviour. 

§3. Variation of -Mat Under a Chiral Invariant Scalar Field 

In order to avoid large violations of chiral symmetry, it is most natural to sup- 
pose that in models such as QMC we are in a representation where ipNipN is chiral 
invariant. For example, if one views the scalar, intermediate range NN interaction as 
the result of two pion exchange this is exactly what one would find. Within such a 
framework the mass of the pion is protected by chiral symmetry, 35 ) even in medium, 
and one expects only very small deviations of the in-medium pion mass from its free 
value - at least for densities typical of normal nuclei. This theoretical expectation is 
certainly supported by recent pionic atom data. 34 ) 

In order to draw conclusions about the scalar polar izability of the nucleon the 
issue to be addressed is then simply: how does the mass of the nucleon change 
when the quark mass varies, m q — > m* = m q — g^cr, with the pion mass kept fixed? 
The extensive studies of the nucleon mass in lattice QCD and especially the success 
in relating full and quenched lattice QCD simulations, 36 -* give us confidence that 
we understand the pion loop contributions. In particular, it was found in Ref. 36 ^ 
that the coefficients were the same, within fitting errors, in full and quenched 
QCD when we used a dipole regulator with mass 0.8 GeV - in the case of both 
the N and the A. We note that the separation of the pion loop contributions is 
model dependent, unlike the extraction of the renormalized coefficients, q, where 
we were careful to establish the model independence at the 1% level. However, the 
possibility of establishing a connection between QCD and fundamental properties 
of nuclear systems, such as saturation, is so important that it seems worthwhile to 
tolerate a degree of model dependence at this stage. 

The free nucleon mass expansion in free space: 

Mn = ao + a^m 2 + a^m\ + self-energy (m n , A) (3-1) 

then becomes 

= ao + a2m* 2 + a^m* 4 + self-energy (m n , A) (3-2) 

in-medium. As noted above, we use a dipole regulator of mass 0.8 GeV for the pion 
loops, which do not change in-medium because the pion mass is protected by chiral 
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symmetry. To lowest order the GMOR relation gives m* 2 oc m*. However, since 
the scalar potential is not necessarily small, we need to keep the next order term 
in order to extract the dependence of the nucleon mass in-medium on the effective 
quark mass. 

In xPT the next term is non-analytic in m q . 37 ^ However, over the relevant mass 
range, fits to lattice data, as well as fits to Schwinger-Dyson phenomenology 38 )' 39 ^ 
and a best fit to the x?T expansion, all suggest that 

ml 4mq + 20m 2 , (3-3) 

with all masses in GeV. 

If we now replace the quark mass, m q , by m* = m q — g$a, to second order in 
a one finds M* N = M N - {Aa 2 gl)a + (20a 2 + 16a 4 )(g^) 2 a 2 , which can be written as 
M^j — Mn ~ <7<t(1 — 9<jv)v-i if we use the numerical values of 02,4 from the analysis 
reported in Sec. 2.2 (without the tadpole). That is, g aa « — g 2 . and the reduction of 
the effective cr-nucleon coupling at nuclear matter density, where g a a £ (150,200) 
MeV is around 15-20 % . This is in remarkable agreement with the upper end of the 
values found phenomenologically in the QMC model. 

If instead of the values found in Sec. 2.2 one used the more recent fit of Table I 
(including the tadpole term), the coefficient of g 2 . in the expression for g aa would 
be -0.6 instead of -1. This is in remarkable agreement with the value given by 
the bag model (g aa = —O.llRg 2 ., with R = lfm = 5GeV _1 ). In terms of the 
phcnomenological role of g aa in nuclear saturation, the difference between these two 
estimates is unimportant. 

To summarize, the FRR expansion of the nucleon mass in terms of quark mass, 
with parameters determined by fits to the best two-flavor, lattice QCD simulations, 
leads to a scalar polarizability of the sign and magnitude found in the QMC model. 
Amongst the natural consequences of this result, as we explained in the Introduction, 
is that nuclear matter will naturally saturate as the density increases. Thus in a very 
meaningful sense the present result connects the saturation of nuclear matter to the 
structure of the nucleon itself within QCD. 

§4. Conclusion 

The remarkable progress in resolving the problem of chiral extrapolation of lat- 
tice QCD data, in the case (such as the nucleon mass) where the chiral coefficients 
are known accurately, gives one confidence that the pion loop contributions are under 
control. In the case of the nucleon one can use this control to estimate the effect of 
applying a chiral invariant scalar field to the nucleon - i.e. to estimate the so-called 
"scalar polarizability" of the nucleon. We stress that, unlike the determination of 
low energy coefficients, which is demonstrably model independent, for the present 
the determination of the scalar polarizability requires a choice of regulator. We chose 
a dipole of mass 0.8 GeV because of the connection between full and quenched QCD 
that has been established in this case. The resulting value is in excellent agreement 
with the range found in the quark-meson coupling model (QMC), which has provided 
an impressive description of many phenomena in nuclear physics. Of most imme- 
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diate significance is the fact that such a scalar polarizability very naturally yields 
nuclear saturation within a relativistic mean field treatment of nuclear matter. Thus 
in a very real sense the results presented here provide a direct connection between 
our growing power to compute hadron properties from QCD itself and fundamental 
properties of atomic nuclei. This is especially interesting in the light of recent ex- 
perimental evidence 43 ) -45 ) that the structure of a bound "nucleon" differs from that 
observed when it is free. 40 ) -42 ) These are important problems that will be more and 
more a focus for our field in the coming years. 
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